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A scenario is presented, based on renormalization group 
(linear perturbation) ideas, which can explain the self- 
similarity and scaling observed in a numerical study of grav- 
itational collapse of radiation fluid. In particular, it is shown 
that the critical exponent /3 and the largest Lyapunov expo- 
nent ReK of the perturbation is related by (5 = (Reft) -1 . We 
find the relevant perturbation mode numerically, and obtain 
a fairly accurate value of the critical exponent /3 ~ 0.3558019, 
also in agreement with that obtained in numerical simulation. 

04.40.-b, 04.70.Bw, 64.60.Ak 



1. Introduction. Gravitational collapse with forma- 
tion of black holes is one of the main problems of classi- 
cal general relativity. Recently, Choptuik 10 discovered a 
"critical behaviour" of the gravitational collapse of mass- 
less scalar field by a numerical study. His result can be 
summarized as follows: Let the initial distribution of the 
scalar field be parameterized smoothly by a parameter 
p, such that the solutions with the initial data p > p* 
contain a black hole while those with p < p* do not. For 
several one-parameter families investigated, near-critical 
solutions (p w p*) satisfy the following: (i) the critical 
solution (i.e. p = p*) is universal in the sense that it ap- 
proaches the identical spacctimc for all families, (ii) the 
critical solution has a discrete self-similarity, and (iii) for 
supercritical solutions (p > p*) the black hole mass satis- 
fies Mbh oc {p — p*) 13 and the critical exponent (3, which 
is about 0.37, is universal for all families. Abrahams and 
Evans || found similar phenomena in axisymmetric col- 
lapse of gravitational wave. 

Evans and Coleman found similar phenomena with 
(3 ~ 0.36 in spherically symmetric collapse of radiation 
fluid, in which case the self-similarity is not discrete but 
continuous. Employing a self-similar ansatz they also 
found a numerical solution which fits the inner region of 
the near-critical solutions very well, and suggested that 
linear analysis around the solution will be useful. 
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In this letter, we present a scenario that explains the 
critical behaviour observed in the radiation fluid collapse. 
We directly relate eigenvalues (Lyapunov exponents) of 
perturbations of the self-similar solution to the critical 
exponent /3, using an argument of renormalization group 
transformation. The formulation is general, and could be 
applied to other models with approximate self-similarity. 
We find the eigenvalues of the perturbation by numeri- 
cal analysis, and find that the value of the exponent (3 
predicted from our analysis matches very well with that 
observed in S. 

We present our picture in sec. 2. After reviewing the 
equations of motion in sec. 3, and the self-similar solu- 
tion in sec. 4, we confirm our picture to an extent in 
sec. 5 by numerical study. Sec. 6 is for conclusions and 
discussions. 

2. A Scenario based on renormalization group ideas. 
We give a formalism for linear perturbations around a 
self-similar spacetime from the point of view of renormal- 
ization group, which proved to be extremely successful in 
the study of critical phenomena [Q. We have benefitted 
from the formulation of g . The argument is general but 
the notation is so chosen as to pave the shortest path to 
our analysis on radiation fluid. We introduce the "scal- 
ing variable" s and the "spatial coordinate" x, which 
are related to the time t and the radial coordinate r by 
s = -ln(-t), x = In(-r/f). 

2a. Renormalization Group and Linear Perturbations. 
Let h — (h\, fi2, h m ) be functions of s and x which 
satisfy a partial differential equation 



dh dh 



o. 



(i) 



Suppose that the PDE is invariant under the "scaling 
transformation" (translation in s) with sq € R 



h(s, x) h- > h(s + sq, x). 



(2) 



A renormalization group transformation (RGT) R SQ is a 
transformation on the space of functions of x, 



where 



1l S0 : H > 



H(x) 
H^ s °\x) 



h(0,x), 
h(so, x). 



(3) 



(4) 
(5) 
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Namely, one obtains by evolving the initial data H 



at s = by the PDE to s 



1Z S0 forms a semigroup 



So- "-so 

with parameter sq, and we denote its generator by D1Z, 
i.e. DTZ = lim (TZ So — 1)/sq. In this context, a self- 

similar solution /i s 



(s,x) — H ss (x) can be considered as 
a fixed point of 1Z S0 for any sq, and is characterized by 
K S0 H SS = H ss or Sftff,, = 0. 

The tangent map of 1Z S at a fixed point H ss is defined 
as a transformation on functions of x: 



T S0 F= lim 

£ — 



Tl 30 (H ss + eF) - H s: 



(6) 



An eigenmode F(x) of DT = 
tion which satisfies (k € C) 



lim (T So — 1)/sq is a func- 



(7) 



These modes determine the flow of the RGT near the 
fixed point. A mode with ReK > 0, a relevant mode, is 
a flow diverging from H ss , and one with Ken < 0, an 
irrelevant mode, is a flow converging to it. A Re/t = 
mode is called marginal. 

2b. The critical solution. In this and the next subsec- 
tions, we present a scenario which explains the observed 
critical behaviour of radiation fluid, assuming that there 
is a unique relevant mode with eigenvalue n (and, for sim- 
plicity, no marginal mode) around the fixed point H ss . 
This assumption is confirmed to some extent in sections 
3-5. 

The assumption implies that the RG flow around the 
fixed point is shrinking, except for the direction of the 
relevant mode (Fig. Q). There will be a "critical surface" 
or a "stable manifold" S of the fixed point H BS , of codi- 
mension one, whose points will all be driven towards _ff ss . 
There will be an "unstable manifold" U of dimension one, 
whose points are all driven away from H ss . A one param- 
eter family of initial data / will in general intersect with 
the critical surface, and the intersection H c will be driven 
to H, s under the RGT: 



lim \H^(x) - H ss (x)\ = 0. 



(8) 



So, H c is the initial data with critical parameter p* . The 
existence of a critical solution for an arbitrarily chosen 
family of initial data thus supports the assumption of a 
unique relevant mode ||. 




FIG. 1. Schematic view of RG flows near H ss . 

2c. The critical behaviour. We now consider the fate 
of an initial data H m i t in the one-parameter family, which 
is close to H c (e = p — p*): 

tfinit(z) = H c (x) + eF(x). (9) 

We evolve this data to s = sq (chosen later): it will 
first be driven towards H ss along the critical surface, but 
eventually be driven away along the unstable manifold. 
Using linear perturbations we have 0]: 



eF) 



~H^ S0) + tf Sa F + O{e 2 ). (10) 

In the second term, only the relevant mode survives: 

f So F = exp( S oST)F~e KSo F rcl , (11) 

where F Te \ is the component of the relevant mode in F. 
Due to (pf), we finally have (for large sq and x < sq) 



HtJ(x) 



, njjl , -H ss (x)+ee KS0 F rcl (x). (12) 

Now we choose so so that the first and second terms 
in ( p"2| ) become comparable, i.e. 

ee (Re«)a =0(1). (13) 

Now that the second term is of 0(1), the data 
differs from H ss so much that one can tell the fate of this 
data depending on the sign of e; and if a black hole is 
formed, the radius of its apparent horizon, and thus its 
mass, will be O(l) measured in x. 

Finally, we translate the above result back into our 
original coordinate (t, r). The relation r = e x ~ s implies 
that the radius of the apparent horizon, which is 0(1) 
measured in x, is in fact O(e~ so ) measured in r. So we 
have from JT3] ) 

M BH ~ O{e~ so ) ~ 0(e 1/(RcK) ). (14) 

Therefore the critical exponent is given exactly by 

1 



Re k 



(15) 
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3. Equations of Motion. The line element of any 
spherically symmetric spacetime is written as 



ds 2 = -a 2 {t,r)dt 2 + a 2 {t,r)dr 2 + r 2 (dd 2 + sin 2 



(16) 



The only coordinate transformation which preserves the 
form of eq. (fLq) is 



(17) 



We assume the matter content is perfect fluid having 
energy-momentum tensor T a b — pu a Ub + p(g a b + u a Ub), 
where u a is a unit timelike vector and p = (7 — l)p. 
We only consider the radiation fluid 7 = 4/3 in this 
letter. The components of u a can be written as ut = 
o(l - V 2 )- 1 ' 2 and u r = aV(l - V 2 )" 1 / 2 , where V 
is the 3-velocity of fluid particles. In terms of vari- 
ables s = — ln(— t),x = ln(— r/t), and introducing N = 
aa~ 1 e~ x , A = a 2 ,(jj = Airr 2 a 2 p, we can write the equa- 
tions of the system in an autonomous form, which makes 
the scale invariance of the system transparent: 



A x „ 2uj ( V 2 
= 1 -A-\ —t 1 1 



A 



1 - V 2 



3 J ' 



N x , 2uo 

— = -2 + A , 

N 3 ' 

u s + (1 + NV)lj iX A{VV s + (N + V)V X } 

w 3(1 - V 2 ) 

NVA X AVN X „ / Ato 

■ \-2N 1 + — — rr =0, 



3A ' 3 ' " V 9(1 -V 2 ) 
AVu), s + (AV + N + -iNV 2 )uj tX 



(18) 



4{(1 + V 2 )V S + {1 + V 2 + 2NV)V, 



1-V 2 



N(l - V 2 )A t , 
A 



+ 4(1 + V 2 )N, X + 2N(1 + W 2 ) = 0. 



4.. The critical solution. One obtains self-similar 
spacetimes by assuming that N and A depend only on 
x: N = N ss (x), A = A ss (x). Conversely, it can be shown 
H that one can express any spherically symmetric self- 
similar spacetimes in that form if one retains the freedom 
of coordinate transformation (|l7j). Then it follows from 
eqs. ( p^|) that uj as and V ss are also functions of x only: 
lu = u> ss (x), V = V ss (x). We fix the coordinate system by 
requiring that the sonic point (see below) be at x = 0. 

We require (i) that the self-similar solution be analytic 
for all x G R, and (ii) as a boundary condition that the 
spacetime and the matter are regular, A = 1 and V = 0, 
at the center (x = — 00). As has been extensively studied 
by Ori and Piran || , the analyticity condition (i) requires 
that the solution be analytic in particular around the 
sonic point, where the velocity of the fluid particle seen 
from the observer on the constant x line is equal to the 
speed of sound 1 / V3- The sonic point is a singular point 



for the ODE's satisfied by self-similar solutions, and con- 
sidering the power series expansion, one can see that the 
solutions are specified by one parameter, say, the value 
of V^ s (0). This, together with the regularity condition at 
the center (ii), restricts T4 S (0) to have only discrete val- 
ues. We employ the Evans-Coleman self-similar solution 
as our H ss in the following. 

5. Perturbation. Perturbation equations (Q) are ob- 
tained by taking the first order variation in eqs. ( |lS| ) 
from the Evans-Coleman solution _ff ss : 



h(s,x) = H ss (x) + e/i var (s,x), 



(19) 



where h represents each of (N,A,u>,V). We require 
■Nvar(S) 0) = to fix the coordinate freedom (|l7|). 

We consider eigenmodes of the form h veir (s,x) = 
h p (x)e KS , with k € C being a constant. Substituting 
this form into (|l9[), and then into ( |l8| ) yields linear, ho- 
mogeneous first order ODE's for (N p , A p , w p , Vp). 

As in the case of self-similar solutions, we require (i) 
that the perturbations are analytic for all icR, and (ii) 
that the perturbed spacetimes are regular at the center 
(A p , V p are finite at x = —00). The sonic point becomes 
a regular singular point for the perturbations. It is not 
hard to see that apart from the overall multiplicative 
factor, the perturbation solutions which satisfy the ana- 
lyticity condition (i) are specified by one free parameter 
k. This, together with the regularity condition (ii) at the 
center, allows only discrete values for k in general. 

Fig. ^ shows the profile of the largest relevant eigen- 
mode obtained numerically. It has the eigenvalue k ~ 
2.81055255; which corresponds to the exponent value 
(3 ~ 0.35580192, according to our scenario of section 2. 
This is in good agreement with the value of 0]. 




FIG. 2. Profile of the eigenmode with the largest eigen- 
value. Curves represent A p (solid line), N BB ~ 1 N P (dashed), 
u!ss~ 1 i^p (dotted), and V p (dot-dashed). 

To further confirm our scenario, we have checked that 
there are no other relevant (or marginal) eigenmodes in 
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the range < Ken < 15, |Iitlk| < 14. We remark that 
there is an unphysical "gauge mode" at k ~ 0.35699 (in 
our gauge) , which emerges from a coordinate transforma- 
tion applied to the self-similar solution. Due to the com- 
plicated structure of the equations of motion, we have 
not found a beautiful argument which can restrict possi- 
ble eigenvalues (like that of |T3] and references therein). 

6. Conclusions and discussions. In this letter we 
have presented a scenario based on the renormalization 
group (linear perturbation) ideas, by which the critical 
behaviour in the black hole formation in the radiation 
fluid collapse is well understood. In particular, we have 
shown that the critical exponent (3 is equal to the in- 
verse of the largest Lyapunov exponent ReK. We have 
performed a partial confirmation of the picture, and mod- 
ulo some assumptions about distributions of eigenmodes, 
have found an accurate value (3 ~ 0.35580192, which is 
also close to that reported in Q . 

Complete confirmation of the scenario requires further 
study of local and global structure of RG flows around 
the fixed point H ss . To establish the local picture, i.e. 
the contraction property of the RGT on the cospace of 
our relevant mode in some neighbourhood of H ss , one 
could try to prove that the eigenmodes form a complete 
set of solutions and that all modes except our relevant 
one are in fact irrelevant. Establishing the global picture, 
which corresponds to proving that the global critical sur- 
face exists and the RG flow around it is as depicted in 
Fig. 1, will pose another challenging problem which will 
be beyond the scope of linear perturbations. Complete 
knowledge of the RG flow could expand the horizon of 
our understanding the gravitational collapse. 

It is expected that the critical behaviour in scalar field 
collapse, where the self-similarity is discrete, can be un- 
derstood in a manner similar to the analysis in this letter. 
One can consider the critical solution as a fixed point of 
the RGT lZs with a suitably chosen so, and can per- 
form linear perturbation analysis around it. A work in 
this direction is now in progress. 

Our intuition on gravitational collapse still seems to be 
heavily based on few exact solutions, especially the lim- 
iting case of pressureless matter. The critical behaviour 
may provide a different limiting case that the final mass 
is small compared to the initial mass for more realistic 
and wider range of matter contents. It will be of great 
help to settle the problems in gravitational collapse such 
as cosmic no hair conjecture. 
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